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Abstract 



O , We numerically obtain a class of soliton solutions for Einstein gravity in (n + 1) dimensions coupled 

to massive abelian gauge fields and with a negative cosmological constant with Lifshitz asymptotic 
behaviour. We find that for all n > 3, a discrete set of magic values for the charge density at the 
origin (guaranteeing an asymptotically Lifshitz geometry) exists when the critical exponent associated 
with the Lifshitz scaling is z — 2; moreover, in all cases, a single magic value is obtained for essentially 
every 1 < z < 2, yet none when z > 2 sufficiently. 

<D ; 1 Introduction 

Since its proposal by Maldacena [I], the AdS/CFT correspondence has proven to be of appreciably broad 
£Nj i theoretical utility, providing new lines of research into both quantum gravity and quantum chromodynamics. 

It conjectures the existence of a holographic duality between strongly interacting field theories and weakly 
coupled gravitational dynamics in an asymptotically AdS bulk spacetime of one dimension greater. 

This idea has been extended in recent years beyond high energy physics to describe strongly coupled 
jy^ . systems in condensed matter physics. In particular, it has enjoyed useful applicability to theories that model 

quantum critical behaviour [3J HI [S] characterized by Lifshitz scaling - that is, a scaling transformation of 
the form 

t -> \ z t, r A -1 r, x -> Ax, (1) 



where z > 1 is a dynamical critical exponent representing the degree of anisotropy between space and time. 
For instance, when 2 = 2, the scaling symmetry given by (Q} is associated with a (2 + l)-dimensional field 
theory of strongly correlated electron systems. 

Such theories are conjectured [6] to be holographically dual to gravitational theories whose solutions are 
asymptotic to the so-called Lifshitz metric, 

ds 2 = I 2 (^-r 2z dt 2 + ^ + r 2 d^j , (2) 

where the coordinates (i,r, x 8 ) are dimensionless and the only length scale in the geometry is I. Metrics 
asymptotic to ([5]) can be generated as solutions to the equations of motion that follow from the action 

1 = lib / dn+lx ^ ( R - 2A - \ F ^ V \ H ^ U - ^ B » B ") . ( 3 ) 

where A is the cosmological constant, F^ v = d^A^ with A^ representing the Maxwell gauge field, and 
Hfj, v = d^BjA is the field strength of the Proca field with mass m 2 = C. 
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This dual theory, referred to as Lifshitz gravity, is known to describe neutral and charged black holes 
[7] [8] [9] [10] whose metrics arc asymptotic to the metric ([2]). It is possible to replace the Proca field with 
higher-order curvature terms |11] | 12 | [T~3"] and attain black hole metrics with the same asymptotic structure. 

Here we explore the possibility of obtaining soliton solutions to the field equations that follow from the 
action © when the Maxwell field is zero. Originally referred to as Lifshitz stars [7], these objects have 
non-singular spacctime geometries with the same asymptotics as their black hole counterparts. We prefer 
to call them solitons since, unlike stars, there is no sharp boundary between a vacuum and non-vacuum 
region. While black hole solutions in (n + 1) dimensions [TU] and soliton solutions in (2 + 1) [T3] and (3 + 1) 
dimensions [7J have already been found, solitons in higher dimensions have not been investigated thus far. 

Here our central aim is to obtain soliton solutions for general (n, z), and to discuss some of their conse- 
quences. For any dimension n, we find that a single soliton solution to the field equations exists for z < 2, 
whereas for z > 2 we are unable to obtain any solutions. For z = 2 we find a discrete set of soliton solutions 
in any dimensionality. Our results are numerical, and so we are able to obtain additional soliton solutions 
for \z — 2 1 < e for sufficiently small e , where e decreases with increasing dimensionality. 

The rest of this paper is organized as follows. In Section 2, we give the field equations for Lifshitz gravity 
and then reduce them to a system of four first-order ODEs. In Section 3, we numerically solve the equations 
of motion for n = 3, 4, 5, 6 and list the corresponding magic values for the central charge density in each case. 
Section 4 concludes the paper, and the Appendix gives the small-radius scries coefficients of the metric and 
gauge functions (needed to determine a set of initial conditions for Section 3). 

2 The (n + l)-Dimensional Field Equations 

Given the action © with = 0, the variational principle yields the field equations [TU]: 



where the energy-momentum tensor of the gauge fields is 



(4) 
(5) 
(G) 



T — — - 



1 



-H pa H f " J g pu — H^Hpv + C ( -B p B p g fll , — B fl B l 



(7) 



The general (n + l)-dimensional metric preserving the basic symmetries (JJ can be written as 



ds 2 



e 2 ( -r 2z f (r) dt 2 + g2 {r \ dr2 + r 2 dn 2 , 



(8) 



where 



n—1 i—l 

^ + EII sin2 ^ e > 2 ' 

i=2 3 = l 

(n-1 i-1 
d6l 2 + II Sin2 6 i d0 i 

i=3 3=2 

n-1 



(9) 



k = 



is the metric of an (n — l)-dimcnsional hypcrsurfacc with constant curvature (n — l)(n — 2)k. 
The Proca field is assumed to be 



B t = q£r z f (r) j (r) , H tr = qtzr^g (r) h (r) / (r) , 



(10) 



with all other components either vanishing or given by antisymmetrization. The asymptotic conditions / (r) , 
g (r), h (r), j (r) — > 1 as r — > oo (required to ensure that any solutions obtained are asymptotic to ([2])) impose 
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the following constraints: 



c= (nzi)£ ?2 = ^1), A 



(z - l) 2 + n (z - 2) + n 2 



(11) 



£ 2 ' y z ' 2£ 2 

It can be shown [TU] that the above reduce the field equations (d))-® to a system of four first-order ODEs 



df 



dh 
dr 



f 



—^-(2 \{n - 1) (z - 1) j 2 - z - 1) h 2 + (z - l) 2 + n (z - 2) + n 2 
1 H I L 



rfr 4(n-l)i 

+ 2 (n - 1) [(n - 2) H 2 g 2 - (n + 2z - 2) r 2 ] |, 



4 (n ^ 1)r2 {2[(n-l)^-l)j 2 +^(^-l)fe 2 -(^-l) 2 -n(^-2)-n s 
- 2 (71 - 1) [0 - 2) fcJ?V - nr 2 ] }, 

=_ 

rf?' 4 (n — 1) 5 

+ 2 (71 - 1) [0 - 2) W 2 g 2 - (n - 2) r 2 ] j + z.gfr, 
(n - 1) (j.g - fc) . 



— J 2 \(n -l)(z- 1) j 2 - z (z - 1) h 2 + (z - l) 2 + n (z - 2) + n 2 
1 r z L L 



2 2 



2 2 
r g 



2 2 
r g 



(12) 

(13) 

(14) 
(15) 



which, in general, cannot be solved analytically. 



3 Numerical Solutions for (n + l)-Dimensional Solitons 

Series solutions for the field equations at large r have been previously obtained [7] [S] [IT] . 

Before we can numerically obtain soliton solutions to the equations of motion (|12[) - p5|) . we require a set 
of initial conditions i.e. the values of /(e), g (e), j (e) and h (e) for some 0<e<l. For soliton solutions we 
demand that the metric be regular and that the Proca field have vanishing field strength at the origin. To 
this end, consider the following series expansions for small r: 



fir) 



1 OC OO OG OO 

-E/f r2P ' 9(r) = rJ29 P r 2p , J (r) = ]>>r 2 ?, ft(r)=r^V 



.2;- 



(16) 



p=0 



p=0 



p=0 



0=1 



It is possible to write all of the coefficients {f p , g p ,j p , h p }°^_ simply in terms of /o and jo (which is the Proca 
charge density at the center of the soliton). In the Appendix, we give the full expressions for the first three 
terms in the series for each function in (fT6|) . Wc find that only for k = 1 are the series solutions finite and 
real, so we shall henceforth set k = 1. 

We can now numerically solve the system (| 12 1) - ([15)) using the shooting method, taking as initial conditions 
(fl~7)l - ([20|) (given in the Appendix) truncated after two terms, evaluated at e = 10~ 6 . Furthermore, we rescale 
all quantities in units of £, effectively setting 1=1. 

As previously noted, (fT7 |) - ([20| are determined solely by two parameters, fo and jo- The value of the 
former can be easily assigned. Since fo appears only as an overall factor in (|T7|) . the system (fT2j) - (fT5|) can 
be solved by setting fo = 1 as an initial condition. If g (r), h (r), j (r) — > 1 and / (r) — >• c 7^ 1 as r — s- 00, wc 
simply change the initial value of / (r),to fo = 1/c and so obtain the desired asymptotic behaviour. 

The determination of the latter, however, is not quite so trivial. We find in general that, as for the 
(3 + l)-dimcnsional case [7], the conditions / (r), g (r), h (r), j (r) — > 1 as r — > 00 can be satisfied only for 
certain discrete values of jo, known as 'magic' values [7]. These correspond to the intercepts of the function 
jo7o; here, 70 := j (r£) — 1, where j (r) is the numerical solution to (fT4"l) dependent upon our choice of jo, and 
Tl is picked to be very large. We furthermore expect these intercepts to coincide with those of the functional 



* Although these intercepts are, of course, the same as just those of oq, /?o and 70, multiplication by jo makes our 
plots more readable. 
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joao and jo^o', here, ao := g (rj,) — 1 and /3o := h (r^) — 1, where g (r) and h (r) are the numerical solutions 
to (fT3|) and (fT5|) respectively, again depending upon the value of jo- For z — 2 in (3 + l)-dimensions these 
intercepts correspond to the removal of a zero mode in the large-r linearized field equations. 

In all of our work, we use = 10 5 unless otherwise stated. Thus, plotting jo7o ( a s a function of jo) for 
different n, and different z for each n, will suffice to give us the magic values, if any exist. 

We will consider three cases: 1 < z < 2, z — 2 and z > 2. 

3.1 Case I: 1 < z < 2 

We find that for all 1 < z < 2 and any n > 3, a single magic value of jo exists. Concordantly, here we 
essentially always observe a single intercept of the function jo7o- The only possible exception is for \z — 2| 
small (generally by less than 0.1) in which case we found numerically that more than one magic value may 
exist. Such scenarios are qualitatively comparable to those where z is exactly 2. Since we discuss this case 
at length in the next subsection we shall not elaborate on the small \z — 2| cases any further here. 

Figure 1 shows this for n = 3 and n = 4, with the situation being qualitatively very similar for all other 
n. Moreover, in Figure 2 we give the magic values that we have computed for n = 3, 4, 5, 6 as a function of 
1 < z < 2. We see that (in any dimension) the magic value increases with increasing z. 




z=1.10 z=1.35 - z=1.60 z=1.85 z=1.10 z=1.35 - — z=1.60 z=1.85 



Figure 1: Plots of jo7o as a function of jo for different values of 1 < z < 2. (The intercepts are the magic 
values corresponding to each case.) Left: n = 3. Right: n = 4. 
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Figure 2: Magic values of jo as a function of 1 < z < 2 for n = 3,4, 5, 6. This plot was obtained by choosing, 
for each n, a discrete set of nine values of z equally spaced along the given interval, numerically computing 
the magic values for each (i.e. the intercepts corresponding to the kinds of graphs depicted in Figure 1), and 
then polynomially interpolating between them. 

3.2 Case II: z = 2 

In this case, a set of magic values exist for any n for all cases we numerically investigated. Figure 3a depicts 
the functions joao, jo/?o and jo7o f° r n = 4, 5 while jo7o for n = 6 is plotted Figure 3b. We find that the 
larger the value of n, the more rapidly oscillating (and hence numerically unstable) these functions become 
- and also, the more difficult we find it to obtain convergence of their intercepts with satisfactory accuracy. 

Figure 4 lists the numerical values thereof for n = 3, 4, 5, 6, i.e. the magic values in each dimension. For 
the plot of jo7o corresponding to n = 3, see [7J; note that in this case, we recover the same magic values as 
those obtained therein. 

Furthermore, we plot the metric and gauge functions corresponding to the lowest magic value for n = 
4, 5, 6 in Figure 5a. Accordingly, all of these are observed to converge to 1 as r — s- oo and, as generally 
expected for solitons, we see that g (r) and h (r) vanish as r — ► 0. In addition, the function r z f (r), again for 
different n, is plotted in Figure 5b. 




Figure 3a: Plots of jo^o, joA) an d jo7o a s functions of jo for z = 2. Left: n = 4. As expected, all three 
functions are observed to have the same intercepts. Right: n = 5 with tl = 10 6 . In this case, convergence of 
the three functions to the same intercepts cannot be obtained numerically quite as accurately. In particular, 
the second intercept (i.e. magic value) is 20.3 for jocto and jo A), but 22.0 for jo7o; the rest, however, are 
found to be the same. 
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jO Y 



1. x 10 < l 



5. x 10 - 7 - 



-5. x 10 7 - 



-1. x 10 (> - 




Figure 3b: Plot of jo7o when n = 6. Numerical convergence of the intercepts of the above plot with those 
of joao and jo A) is not very accurate in this case (becoming even worse for n > 6), and so here we base the 
list of magic values given in Figure 4 solely on the former. 



n 


Magic values for j when z = 2 


3 


3.59, 21.8, 1.34 x 10 2 , ... 


4 


3.80, 18.0, 91.2, ... 


5 


4.60, 20.3, 1.01 x 10 2 , ... 


6 


7.20, 9.74, 13.5, ... 



Figure 4: The three lowest magic values for z = 2 and n = 3, 4, 5, 6. 




■f(n=4) 
f(n=6) 



■g(n=4) 
g(n=6) 



f(n=5) 



-g(n=5) 




■ h (n=4) 
h (n=6) 



-j(n=4) 
J (n=6) 



h (n=5) 



-j(n=5) 



Figure 5a: The metric and gauge functions corresponding to the lowest magic value for n = 4, 5, 6. Left: 
The metric functions / (r) and g (r). Right: The gauge functions h (r) and j (r). 
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r 

(r A z)*f(n=4) (r A z)*f(n=5) (r A z)*f(n=6) 



Figure 5b: The function r z f (r) for n = 4, 5, 6, illustrating that we indeed have soliton solutions. 
3.3 Case III: z > 2 

While (multiple) intercepts of the function jo7o can still be obtained if z is made only very slightly greater 
than 2, none are seen to occur, in any dimension, once z becomes appreciably greater. In other words, we 
found no more magic values for z > 2. In particular, none exist for the zero modes (i.e. for z = n — 1) of 
any n > 3. In Figures 6a, 6b and 6c we plot jo7o for, respectively, n = 4, 5, 6 and various z > 2. We find 
that magic values can still be obtained (i.e. intercepts of these graphs exist) only if, for instance, we have 
approximately z < 2.021 when n = 4, z < 2.005 when n = 5, and z < 2.002 when n — 6. The analogous 
plots of joQ^o and joA) are found to be qualitatively similar, and so we have not included them here. 

J -0.01 -n 1 1 1 1 , r- 




z=2.0001 z=2.001 — z=2.01 z=2.1 z=3 



Figure 6a: Plots of jo7o as a function of jo for n = 4 and various z > 2. 
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z=2.0001 z=2.001 — z=2.01 z=2.1 z=3 z=4 



Figure 6b: Plots of jo7o as a function of jo for n = 5 and various z > 2. 




z=2.0001 z=2.001 — z=2.01 z=2.1 z=3 z=4 z=5 



Figure 6c: Plots of jo7o as a function of jo for n = 6 and various z > 2. 

4 Conclusion 

We have searched for soliton solutions in asymptotically Lifshitz spacetimes from (3 + 1) [7] to (n + 1) 
dimensions. We have found that such solutions do exist, but with a somewhat surprising consequence: 
namely, solutions exist only when the critical exponent associated with the Lifshitz scaling is (very close to) 
2, or smaller. In particular, 1 < z < 2 and z « 2 (to within at most 0.021, Vn > 4) yield, respectively, a 
single magic value and a discrete set of magic values, in any dimension. But, once z exceeds 2 sufficiently, no 
more magic values - and hence no more solutions for the metric and gauge functions - can be numerically 
found. This means, therefore, that no soliton solutions exist for the zero modes of any n > 3. 

It would be interesting to understand in greater depth why z = 2 is such a special point in parameter 
space in all of the dimensions we investigated. Are such solutions stable, or will they undergo collapse to a 
black hole? The relationship between these solutions and the general (in)stability of asymptotically Lifshitz 
spacetimes [15] would be another interesting subject to investigate. 
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Appendix 



It can easily be verified that 



where 
h 
h 



fa 



2k£ 2 n{n- 1) 
-fo 



f(r) = -(/o + /ir 2 +/ 2 r 4 ) ; 

g(r) = r (g + g x r 2 + g 2 r 4 ) , 

j{r) = jo +jir 2 + j 2 r 4 , 

h(r) = r (h + h x r 2 + h 2 r 4 ) , 



{ [(« - if + z (n - 2) + z 2 ] + [(z - 1) (n - l) 2 ] j 2 } , 

2(z- l) 4 + (z 2 + 2z - 7) (z - l) 2 n + 2 (z 3 - z 2 - 3z + 4) n 2 



(17) 

(18) 
(19) 
(20) 



8£ 4 k 2 (n - iy n 2 (n + 2) 

+ (3z 2 - 4z - 2) n 3 + 2 (z - 1) n 4 + n 5 j + (z - 1) [4 (z 2 - 3z + l) - 2 (z 2 - lOz + 5) n - 4 (z 2 + z - l) n 2 
+ 2 (z - 1) (z - 4) n 3 + 8 (z - 1) n 4 - 2 (z - 1) n 5 |j^ + (z - l) 2 [2 - 3n - 4n 2 + 10n 3 - 6n 4 + n 5 ] j 4 



WW 



91 2k 3 / 2 £ 3 n (n - 1) 
-1 

.92 



8^ 5 fc 5 / 2 (?i-l) (n + 2) 



{ [n 2 + 1 - (2 - z) („ + z)] - [(z - 1) (n - 1)] j 2 } 



6 (z - l) 4 - 3 (z 2 +2z-7)(z- l) 2 n + (-6z 3 + 6z 2 + 18z - 24) ■ 



+ (-9z 2 + 12z + 6) n 3 + 6 (-z + 1) n 4 - 3n 5 + 2 (z - 1) - 2 (3z 2 - 7z + 3) + (z 2 - 12z + 13) 



+ (5z 2 - 9z - 3) n 2 + 9 (z - 1) n 3 + (-2z + 5) n 4 j 2 + (z - l) 2 (n - l) 2 (4?i 2 - 7n - 6) 



.72 



-Jo 



2fc€ 2 n(n- 1) 
jo 



{ [(« - I) 2 + z (3n - n 2 + z - 3)] + [(z - 1) (n - l) 2 ] Jo) , 

2 (3z 2 - 8z + 3) (z - l) 2 + (3z 4 + 4z 3 - 41z 2 + 60z - 2l) 



S£ 4 k 2 (n-l)n 2 (n + 2) 

+ (lOz 3 - 14z 2 - 26z + 24) n 2 + (-4z 3 + 27z 2 - 20z - 6) n 3 + (-8z 2 + 18z - 6) n 4 + (z - 1) (z - 3) n 5 
+ 2 (z - 1) [2 (3z 2 - 8z + 3) - (7z 2 - 31z + 19) n - 2 (z 2 + 5z - 9) n 2 + z (3z - 13) n 3 + 2 (5z - 4) n 4 



2 (2z - 3) n 5 j 2 + (z - l) 2 (n - l) 3 (3n 



6) 
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ho = 



hi 



jo (n- 1) 
-Jo 



"2fc 3 / 2 £ 3 n (n + 2] 



{ [n 2 (2 - z) + 4n (z - 1) + 2z 2 - 5z + 2] + [(z - 1) (n 2 - 3n + 2)] 



Jo 



8 (2z - 1) - 2) (z - if + (8z 4 + 10z 3 - 107z 2 + 154z - 56) 



8^ 5 fc 5 / 2 (n - 1) (n + 4) n 2 (n + 2 
+ (20z 3 - 18z 2 - 72z + 64) n 2 + (-6z 3 + 46z 2 - 36z - 16) n 3 + (-10z 2 + 32z - 16) n 4 + (z - 2) (z - 4) n 5 

+ 4 (z - 1) U(2z - 1)0 - 2) - 2 (3z 2 - 15z + ll) n - (4z 2 + 3z - 16) n 2 + (2z 2 - 13z + 4) n 3 + (7z - 8) n 4 
+ (-z + 2) 7i 5 j j 2 + (z - l) 2 Un - l) 2 (n - 2) (3?i 2 - 4n - 8)1 

satisfy the four ODEs (fT2l) -(fl~5 | . at least up to sixth order in r. In Figure 7, we compare the above scries 
solutions with the numerical results obtained (in Figure 5a) for the lowest magic value when n = 4, and find 
that they are in good agreement up to at least r w 0.1. These plots are found to be qualitatively similar for 
all other n > 4, and hence we refrain from showing them as well. 




1.5 



■ f series (n=4) 
• g series (n=4) 



f numeric (n=4) 
g numeric (n=4) 




- h series (n=4) 

- j series (n=4) 



■ h numeric (n=4) 
" j numeric (n=4) 



Figure 7: The metric and gauge functions corresponding to the lowest magic value for n = 4: numeric 
solutions (Figure 5a) versus series solutions ((fT7 j) -(|2"0" |) . truncated after two terms) for small r. Left: The 
metric functions / (r) and g (r). Right: The gauge functions h (r) and j (r). 
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